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We report Monte Carlo studies of a two-dimensional soft colloidal crystal confined in a strip ge- 
ometry by parallel walls. The wall-particle interaction has corrugations along the length of the strip. 
Compressing the crystal by decreasing the distance between the walls induces a structural transi- 
tion characterized by the sudden appearance of a one-dimensional array of extended defects each 
of which span several lattice parameters, a "soliton staircase" . We obtain the effective interaction 
between these defects. A Lindemann criterion shows that the reduction of dimensionality causes a 
finite periodic chain of these defects to readily melt as the temperature is raised. We discuss possible 
experimental realizations and speculate on potential applications. 



PACS numbers: 74.25. Qt,61.43.Sj,83.80.Hj,05.65.+b 

There are many examples of condensed matter systems 
where extended defects in some order parameter field be- 
have as effective "particles" which themselves undergo 
order-disorder transitions with important consequences 
for the properties of the original system. One can easily 
recall many examples such as charge [1 or spin [2] den- 
sity waves, vortex matter[3 j , Skyrmions[4 in fractional 
quantum Hall systems, domain walls in commensurate- 
incommensurate phases etc. In many of these exam- 
ples, the typical size of these defects is much larger than 
the smallest relevant microscopic length scale. Investi- 
gation of the properties of such defect lattices requires 
knowledge of the effective interactions between defects 
which are usually difficult to measure directly in exper- 
iments. They are also difficult to obtain from computer 
simulations because of the large difference in length scales 
involved and can usually be computed only within a mean 
field approach and in the highly dilute limit [6]. In this 
Rapid Communication we describe a simple example in- 
volving extended defects in a colloidal solid pj [5] where, 
on the other hand, such effective interactions may be ob- 
tained to great accuracy using a relatively small system 
with appropriate use of finite size techniques. 

In a recent work [9 , we have shown that one can 
produce novel defect states in a colloidal crystal con- 
fined within a narrow quasi one dimensional strip by de- 
forming it in a suitable way. We perform Monte Carlo 
simulations [10 of a simple model solid with particles in- 
teracting with a potentialHU [12] V(r) = e(^) 12 at a 
distance r. The solid consisting of n y x n x unit cells with 
lattice parameter a is confined in a channel of length 
L y = riya and width D = (n x — A)av // 3/2, where A is 
a "misfit" parameter (Fig jl] (inset)). Periodic boundary 
conditions are assumed in the x direction whereas, in 



the y direction, the crystalline strip is confined by two 
fixed walls composed of two rows of immobile particles. 
When A = 0, one obtains a triangular crystalline solid 
between the walls at zero tensile stress a. With increas- 
ing misfit A (i.e. tensile strain) a increases up to some 
critical value, where a transition occurs that reduces n x 
by one. At constant density, the n y extra particles of 
the row that disappears are added to the n x — 3 inner 
rows of the strip; the resulting average lattice spacing 
a' = a(n x — 2)/{n x — 3) is incommensurate with the ef- 
fective periodic potential due to the rows of fixed parti- 
cles. This leads to the formation of a "soliton staircase" 
[9] along the length of the walls, (accompanied by a pat- 
tern of standing strain waves in the crystal) according to 
the Frenkel-Kontorova mechanism[13 . The number of 
solitons produced is given by n s = n y /(n x — 3) [9] since 
each soliton contais just a single excess particle. Here 
we extend our study and investigate the structural and 
mechanical properties of the system, and show that the 
soliton superstructure in confined crystals behaves as a 
one dimensional system of extended "particle" -like exci- 
tations which interact among themselves via an "effec- 
tive" harmonic potential. We show how to extract the 
harmonic "spring constant" of this effective lattice and 
study the gradual melting of the soliton lattice into a soli- 
ton fluid caused by raising the temperature. We expect 
our calculations to be of direct relevance to experiments 
on confined colloidal crystals [7\ 151 fT4]. 

In order to extract properties of the defect system, we 
need to obtain the size and position of the individual de- 
fects. We describe below the two independent techniques 
that we have used for this purpose. We shall also show 
that the results obtained by these two methods agree 
with each other making us confident of our conclusions. 
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FIG. 1: Internal stress a — a xx — cr yy (in Lennard- Jones 
units) in the confined crystalline strip plotted vs. A, for the 
case of a system started with n x = 30, n y = 108 (full sym- 
bols) and a system started with n x = 29, n y — 108 with 
the 108 extra particles distributed among the 27 inner rows, 
(open symbols). Curves are guides to the eye only. Inset: A 
schematic sketch of our geometry: we study a system of size 
D in x-direction and L y in y-direction, apply periodic bound- 
ary condition along the ?/-axis, while the boundary in the 
x-direction is created by two rows of fixed particles (shaded) 
on the ideal positions of a perfect triangular lattice at each 
side. In the fully commensurate case, D = n x aV3/2. 



FIG. 2: Characterization of the 6-membered SP rings: (a) 
the modified 6-membered ring where C is the central particle, 
(b) an SP-6-membered ring containing three antipodal pairs 
(al, a2), (bl, b2), (cl, c2) the shortest paths of equal distance 
between (al) and (a2) are through (bl) and (cl) or through 
(c2) and (b2). (c) The results of two-dimensional ring analy- 
sis on 108 x 30 colloidal crystals. In each case, a defect lattice 
consisting of n s = 4 = 108/(30 — 3) solitons was stabilized 
at temperatures T = 0.7(i), 0.5(H), 0.3(m) and 0.1(iv). Note 
that as the temperature is increased, random thermal fluctua- 
tions produce defects within the bulk of the strip, making the 
identification of the defect lattice more and more ambiguous. 



1. The ring method: Here, we identify the particles 
which belong to a single defect using the topologically 
defined concept of shortest path (SP) ring structures; the 
procedure is discussed at length in Refs. [I5j HU [17], to 
which we refer the interested reader for details. Briefly, 
we classify particles as belonging either to a standard 6- 
membered SP ring or to a modified 6-membered SP ring. 
A particle belongs to an SP ring if the number of bonds 
passed through in moving from one particle of the ring to 
another is the shortest among all possible paths through 
the network of bonds. If, in addition, every particle of 
the ring is also bonded to a single central particle, the 
particles are said to belong to a modified SP ring. By 
definition, particles belonging only to SP rings represent 
regions containing defects, while those belonging to mod- 
ified SP rings represent ideal crystalline arrangement. 
Once the positions of the atoms belonging to defect (soli- 
ton) locations are obtained, one can use standard clus- 
ter counting techniques to obtain the coordinates of the 
atoms associated with each individual defect. We observe 
that (1) the defects are extended structures consisting of 
more than one atom and (2) the number of atoms com- 
prising a defect is more or less fixed. We can then easily 
obtain the center of mass coordinates of each defect (i.e. 
group of pink particles, see Fig. [2^c)). In Fig{2|c), the 
red and pink colors represent the locally ideal and defec- 
tive neighborhoods, respectively, inside the 2D strained 
colloidal crystal showing the stabilization of a soliton lat- 



tice at four different temperatures. Note that, as the 
temperature is increased, the location and size of the de- 
fects become ambiguous due to the presence of random 
thermal fluctuations. 

2. The block method: It is important to make sure that 
the soliton positions that we extract from the simulation, 
are not affected seriously by the method used to analyze 
these configurations. Therefore we have used another 
method, viz. the blocking method, to obtain an inde- 
pendent estimation. Accordingly, we introduce a block 
length L\y — n^a, where a is the lattice constant of the 
undeformed lattice. We choose = 8 — 12 and hence 
L5 such that ^> 1 but L5 still clearly less than the 
expected value of At/ • By moving this coarse-graining 
block along a row in the ^-direction (in steps of a), we can 
count how many particles n actually fall inside a block. 
If we work at low enough temperatures, where the mean- 
square displacement of the particles on the lengths L5 is 
still clearly much less than a 2 , we obtain n = if no soli- 
ton core falls into the block, while we obtain n = + 1 
if a soliton core falls into the block. Calculating then 
the center of mass of a cluster of adjacent blocks with 
n = rib + 1 then yields an alternative estimate for the 
position of a soliton in a system configuration. 

It is now straightforward to obtain both the structure 
factor S 80 i(q) = ^ J2u> ( ex P[i Q (Vl ~ Vi')]) and the prob- 
ability distribution P(Ayo) of Ayo the distance between 
the centre of mass positions yi,yi+i of neighboring soli- 
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FIG. 3: Structure factor of the soliton staircase S so i(q) vs. q 
at temperatures T = 0.1(+), 0.2(*), 0.3 (o), 0.4 (A) and 0.6 
(□). Symbols are the actual simulation data obtained from 
either the ring-analysis (a) or the block analysis (b), lines are 
guides to the eye. Each curve has been shifted by an unit 
distance in the vertical direction for clarity. All data are for 
the system with L = n y a with n y = 108 , and n s = 6 solitons 
present after a transition from n x — 20 to n x — 19. Note that 
both (a) and (b) give comparable results. 



tons 1,1 + 1. The sum runs over the positions of all soli- 
tons in a row, and the result is averaged over all rows over 
both boundaries of the system apart from the statistical 
average over many independent particle configurations). 
While at T — in the limit L y and n s — > oo the soliton 
lattice would cause sharp Bragg peaks, for all T > the 
soliton system is expected to have a liquid-like structure 
factor. In fact, assuming that the interaction between 
neighboring solitons is harmonic, 



s «*(q) = ^J2 e ^( l - l ') A yo]^M-\q 2 ((Ui-u ll ) 2 )) 

x w 

(1) 

Each soliton is described as an effective point particle 
of mass M, position yi and conjugate momentum II/. The 
Hamiltonian for the harmonic chain is given by, 

H sol = \ Y^i/M + MC 2 (y l+1 -yi- Ay ) 2 / (Ay ) 2 } 
i 

(2) 

Where the parameter C plays the role of a sound velocity. 
From Eq. [2] one obtains the correlation function of the 



mean square displacements JJ\ = y\ — I Ay as [T8] 

m - U ) 2 ) = l(Ay ) 2 k B T/(MC 2 ) 

= IS 2 (3) 

Here 5 characterizes the local displacement, S 2 = 
((Ui+i — Ui) 2 ) = ((Ayi) 2 ). The obvious interpretation 
of Eq. [3] is that the relative displacements Ayi at each 
index I of the one-dimensional soliton lattice add up in a 
random- walk- like fashion [T2]. From Eqs. Tj3 1 one may 
derive S so i(q) in one dimension, for n x — > oop~8j [T9] 
exactly. However, our data corresponds to very small 
n s ~ 4 — 6. Nevertheless, following Ref. [19 one can 
evaluate Eq. [I] for finite n s , using S and Ayo as parame- 
ters, 
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- cos(qlAy )exp(--q 2 5 2 l) 



(4) 
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As long as S = ((Ay) 2 ) 1 / 2 <C Ayo, the one-dimensional 
correlation extends over many solitons, and the term 
soliton lattice still is in a sense meaningful; when 5 is 
no longer much smaller than A^o, however, the sys- 
tem rather should be described as a soliton liquid. As 
is well known, the melting of a one-dimensional crystal 
is a continuous transition (£ ~ 1/T). If one nevertheless 
defines [17] an effective melting temperature T^ 1 for one- 
dimensional systems by arbitrarily requiring that their 
Lindemann parameter 5 2 / Ayo < 0.01, one would obtain 
ksT^ 1 = 0.0 IMC 2 as the temperature scale that con- 
trols the melting of the soliton lattice. Though clearly 
the melting of the soliton lattice is far from a sharp ther- 
modynamic phase transition, this already suggests that 
the soliton lattice may melt at rather low temperature, 
far below the melting temperature T m 1.35 of the bulk 
two-dimensional crystal at the chosen density [TT]. 

Fig. [3] shows simulation data for S so i(q) vs. q at var- 
ious temperatures. We note that at T = 0.1 indeed 
the peaks of S so i(q) are already rather sharp, while for 
T > 0.3 the structure factor clearly has the character of 
a fluid. The nature of the curves are well represented by 
the form given in Eq(4]and we may obtain a value for the 
(only) parameter 5 by fitting the data. However, we use 
the more accurate procedure discussed below. 

We now compare the probability distribution P(Ayo) 
of the distance between neighboring solitons in the lat- 
tice from both the ring analysis method and this block- 
ing method. Fig. [4] shows that both methods of ana- 
lyzing the configurations to identify where the solitons 
are, agree almost perfectly with each other, and more- 
over P(Ayo) is nicely described by a Gaussian, P(Ay ) = 
Aexp( — I b (Ay — Ayo) 2 ) where A is a constant ensuring 
normalization, Ayo = 18 = n y /n s = 108/6 has been used 
to be consistent with our choice of L and D (an uncon- 
strained fit gives Ayo =, 17.866 ± 0.023 !) so the only fit 
parameter is b = (Ay 2 ) -1 . From Eq. [2] it is obvious that 
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FIG. 4: Probability distribution P(Ay ) versus Ay , for T = 
0.1(a), 0.2(b), 0.3 (c) and 0.6 (d). Open and closed symbols 
are the results from the block and the 2d ring analysis method 
respectively. The lines show the corresponding gaussian fits. 



the quantity K = MC 2 / (Ayo) 2 plays the role of an elas- 
tic constant such that b = K/ksT. FigJE] therefore plots 
ksT/K vs. T, to test to what extent K is independent 
of temperature. Finally, we note that the Lindemann ra- 
tio 0.01 is reached at an effective melting temperature 
Tm l ~ 0-44. This value is compatible with the direct ob- 
servation of the effective melting temperature obtained 
from the structure factor. 
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FIG. 5: Plot of k B T/K, as extracted from P(Ay ), versus T. 
Symbols have the same meaning as in Fig. 4. The straight 
line is a fit to all data, yielding K — 7.3 ±0.3. 



In our previous studies [9], we briefly reported our dis- 
covery of the soliton staircases and the strain wave pat- 
terns in confined soft two-dimensional colloidal crystals of 
different sizes. In the present paper, we obtain accurate 
effective interactions between the solitions and report a 
gradual melting of the soliton superlattice in analogy to 
that of harmonic chains. We believe that our studies 
would be useful in designing experimental colloidal sys- 



tems where such defect lattices in narrow channels may 
be stabilized. These structures, which have some resem- 
blance to vortex matter in channels, [20] may have inter- 
esting optical and transport properties. Work along these 
lines are in progress. 
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